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ABSTRACT 

Bright single pulses of many radio pulsars show rapid intensity fluctuations (called 
microstructnre) when observed with time resolutions of tens of microseconds. Here, we 
report an analysis of Arecibo 59.5 p,sec-resolution polarimetric observations of 11 P-band 
and 32 L-band pulsars with periods ranging from 150 msec to 3.7 sec. These higher 
frequency observations forms the most reliable basis for detailed microstructnre studies. 
Close inspection of individual pulses reveals that most pulses exhibit quasiperiodicities 
with a well-deflned periodicity timescale (P^). While we And some pulses with deeply 
modulating microstructnre, most pulses show low-amplitude modulations on top of 
broad smooth subpulses features, thereby making it difficult to infer periodicities. We 
have developed a method for such low-amplitude fluctuations wherein a smooth subpnlse 
envelope is subtracted from each de-noised subpulse; the fluctuating portion of each 
subpnlse is then used to estimate via autocorrelation analysis. We find that the 
microstructnre timescale P^ is common across all Stokes parameters of polarized pulsar 
signals. Moreover, no clear signature of curvature radiation in vacuum in highly resolved 
microstructnres was found. Our analysis further shows strong correlation between P^ 
and the pulsar period P. We discuss implications of this result in terms of a coherent 
radiation model wherein radio emission arises due to formation and acceleration of 
electron-positron pairs in an inner vacuum gap over magnetic polar cap, and a subpnlse 
corresponds to a series of non-stationary sparking discharges. We argue that in this 
model, P^ reflects the temporal modulation of non-stationary plasma flow. 


Subject headings: MHD — plasmas — pulsars: general, radiation mechanism: nonther- 
mal 


1. Introduction 

Soon after the discovery of pulsars, Craft et al. (1968) noticed snb-millisecond (be., hundreds 
of microseconds) intensity structures in single pulses of PSR B0950+08. Subsequently, several high- 
time-resolution studies revealed rapid intensity variations superposed on the subpulses of several 
other pulsars. This phenomenon, usually known as pulsar microstructure, is currently thought to be 
associated closely with mechanisms of magnetospheric radio emission. Observationally, microstruc¬ 
tures are found to occur simultaneously over a wide band {e.g., Rickett et al. 1975) and in pulsars 
with rotation periods down to 89 msec (Kramer et al. 2002), providing strong empirical evidence 
in favor of magnetospheric emission. The high brightness temperatures observed in pulsars rule 
out any thermal plasma processes, and indicate that the underlying mechanism must be a coherent 
plasma process. However, a coherent radiation mechanism that provides a physical interpretation 
of microstructnre behavior is yet to be identified. 

Broadly, two appealing scenarios have been considered to explain microstructnre emission, 
commonly known as angular beaming or temporal modulation. The angular beaming scheme is 
a geometrical phenomenon wherein intensity fiuctnations result from an angular radiation pattern 
(mimicking microstructnre width and periodicity) in the transverse direction. As the pulsar rotates, 
an observer then samples this emission structure. On the other hand, the temporal modulation 
scheme entails intensity fiuctnation of propagating waves in the magnetosphere that finally reach 
the observer. 

The existing microstructnre studies fall into two broad categories: Studies with very high time 
resolution (~100 nanosec) which are needed to study narrow microstructures; and those with a 
coarser resolution (>10 fisec) that are useful for exploring the properties of broad microstructures. 
Amongst the former is the well-known detection of ultrafast 2 -nanosecond temporal structures in 
the Crab pulsar (Hankins et al. 2003; see also Jessner et al. 2010). Intensity fiuctnations with time 
scales as small as 2.5 //sec have also been reported in PSR B1133+16 by Bartel & Hankins (1982), 
and Popov et al. ( 2002 ) find micropulses with widths between 2 and 7 //sec in pulsars B0950+08, 
B1133+16 and B1929+10. Most of the other well-known microstrnctnre studies were geared to 
study the broad microstructures. These have shown that subpulses often show qnasiperiodic in¬ 
tensity variations with periods of hundreds of microseconds (e. 5 '., Hankins 1971; Boriakoff 1976; 
Soglasnov et al. 1981; Soglasnov et al. 1983; Cordes et al. 1990). Cordes (1979) and Kramer et al. 
( 2002 ) further show that the measured temporal width (t^) of a micropulse is related linearly to the 
rotation period P approximately as (msec) ~ 10“^P (sec) and (msec) ^ 0.7 x 10“^P (sec) re¬ 
spectively. Polarized microstrnctnre has been investigated for a few pulsars (Ferguson & Seiradakis 
1978; Cordes 1979), and only for the Vela pulsar B0833-45 does some quantitative description exist 
(Kramer et al. 2002). Further observations, analyses, and phenomenological assessments are needed 
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to establish how microstructure is couuected physically with the mechauisms of pulsar emissiou. 

Oue of the difficulties iu usiug the existiug microstructure studies to amalgamate with pulsar 
emissiou theories is that they represeut differeut resolutious, frequeucies of observatiou, aud aualy- 
sis methods. Heuce, iu the remaiuder of our work, our maiu effort is to (1) carry out a systematic, 
siugle-frequeucy, siugle iustrumeut (Arecibo) study of the properties of broad polarized microstruc¬ 
ture; aud ( 2 ) study the physical implicatious of the microstructure characteristics iu the coutext of 
curreutly-favored cohereut radio emissiou theories. 


2. Observations 

High-time-resolutiou siugle-pulse polarimetric observatious were carried out usiug the Arecibo 
telescope for 32 pulsars iu the frequeucy rauge of about 1.2 - 1.6 GHz (L-baud), aud 11 pulsars at 
arouud 327 MHz (P-baud). For these observatious, we used the Observatory’s Mock spectrometers 
iu the “siugle-pixel” polarimetry mode (see http://www.naic. edu/$\sim$astro/mock. shtml). 
These Fourier-trausform-based iustrumeuts are capable of sampliug at ~59.5 //sec rate aud produc- 
iug Stokes spectra iu a hlterbauk format. The details of the observatious are giveu iu Table 1. This 
moderate time resolutiou was choseu both as suitable for iuvestigatiug broad microstructure aud 
maximiziug the uumber of objects that could be so observed. Iu the remaiuder of this paper, we will 
use the word microstructure to imply broad microstructure uuless specihed otherwise. Heretofore, 
ouly some 14 weh-kuowu bright pulsars have beeu subjected to microstructure studies (see Table 2 
of Kramer et al. 2002 aud the Vela Pulsar), out of which 11 are iu the Arecibo decliuatiou rauge. 
Our sample iucluded 10 of these previously studied objects (except PSR B0540+23) at L baud as 
well as 22 others that are uow accessible owiug to both the Gregoriau feed system aud capacious 
Mock back-euds. At P baud, there is au overlap with 6 pulsars; however, for most pulsars at this 
frequeucy (apart from B1929+10), either the dispersiou or the scatteriug timescale is comparable 
with the sampliug rate aud heuce they are uot suitable for microstructure study. Nouetheless, 
giveu the uovelty of this 59.5-//sec polarimetry with exceheut sigual to uoise (S/N) ratio, we do 
display their properties iu the ouliue supplemeutary material. The high-time-resolutiou average 
polarizatiou prohles of 24 pulsars at L baud which were eveutually used for microstructure study 
(see Sec. 3) are showu iu Sec. B, while the rest are available iu the ouliue supplemeutary mate¬ 
rial. The pulsar positiou augles are rotated to iuhuite frequeucy usiug the rotatiou measure values 
giveu iu Table 1 aud, iu priuciple, are the absolute positiou augles. We uote that the absolute 
positiou augles of our observatious at L baud are more robust aud the positiou augles appear to 
be iu good agreemeut with earlier published results (e. 5 -., Morris et al. 1981; Johustou et al. 2005; 
Raukiu 2007). Iu P baud, the baudwidths used for our observatious were uot large euough to 
calculate accurate rotatiou measures, aud heuce au error iu positiou augles of about 1 rad m“^ 
(which cau result due to the iouosphere) cauuot be ruled out 


^ Note, that the P band observations will therefore not be included in the analyses below or discussed therein. 
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The average polarization properties of the pulsars in our study agree with lower-resolution 
observations reported in several other studies (e. 5 -., Everett & Weisberg 2001). None of these av¬ 
erage profiles show any obvious signatures of microstructure, and are consistent with the view 
that, at any given longitude, the averaging process washes out any such structure that can be 
seen in individual pulses. The pulse-phase-averaged percentage linear and circular polarization 
values over bins for signal to noise ratio greater than 3 (where the noise was estimated as off 
pulse rms values for linear and circular polarization respectively) for all our observations are given 
in Table 1. The pulse window or phase region over which the averaging was done corresponds 
to outer edges of the pulse profile at three times the off pulse rms level of the total intensity In 
Table 1, we have not included error estimates on %L and %V. RMS values could be computed 
using standard methods of error propagation [see, e.g., Nontsos et al. (2015)] assuming Gaussian 
distributions. However, rigorous error estimates are best-quantified in terms of valid confidence 
intervals, and we believe that such computation for %L oc L/I and %V oc V/I may not be 
straightforward. While the noise in the measured Stokes parameters I^Q^U^V is Gaussian, the 
distribution of the ratio of two Gaussian quantities (i.e., V/I) is not Gaussian (e.g., the ratio of 
two mean-zero Gaussian quantities has a Cauchy or Lorentz distribution; and the distribution of 
L = y/o^TiP is skewed and non-Ganssian such that the ratio %L is non-Ganssian as well). 
Given the non-Ganssianity of %L and %V, characterizing their errors as an RMS serves no pur¬ 
pose. Stokes parameter measurements are also subject to systematic errors due to cross-coupling 
in the feed system. The crossed-linear P and L-band systems at Arecibo have isolations in voltage 
of some 25 db or better (see Arecibo technical memos in http://www.naic.edu/~astro/aotms/ and 
http://www.naic.edn/~phil/stokes/stokescal.html), or in power of about 0.5%. Our observations 
correct for the instrumental gains and phase using a linearly polarized calibration signal, but no 
correction was made for the cross-coupling. We except this latter to introduce a maximum mutual 
error of 5% in L and V (e.g., AL ~ 0.05E and AV ~ 0.05L). 

Single pulses, however, exhibit highly complex structures, with intensity fluctuations over 
timescales much smaller than the widths of their subpulses. Some obvious conclusions can be drawn 
from a visual inspection of individual pulses. For example, in several pulsars, we And that there 
are occasional bright pulses where the subpulses show a highly modulated ringing structure with 
all the Stokes parameters fluctuating with the same periodicity. Several such examples are shown 
in Fig. 1. However, almost all snbpulses show low-amplitude qnasiperiodic intensity variations 
superposed on a broad smooth envelope. In the subsequent sections, we will develop and apply 
a method of extracting microstructure periodicities for both low-amplitude and highly-modulated 
snbpulses. The microstructure emission is also seen in both the orthogonal polarization emission 
modes (an example is shown in Fig. 2). 
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3. Estimating the microstructure timescale for an individual pulse 


The focus of our analysis is on extracting microstructnre timescales from time-series data for 
individual pulses. Timescales are traditionally identified as minima, maxima, or slope changes in 
the autocorrelation function (ACT) for a pulse. We define the (normalized) ACT at lag r for a 


signal I{t) as 


ACF(t) = 


/ + T)dt 

J P{t)dt 


SO that ACF(O) = 1. The behavior of the ACF needs a closer look to understand what can and 
cannot be inferred about the timescales of a signal. For illustration, consider the toy train of 
Gaussian spikes (Fig. 3) which is characterized by two timescales; namely, the baseline width of 
each spike, and the distance between two successive spikes, be., the period of the spike train. 
The gap g = between two successive spikes is assumed to be non-negative. When < g, 

it can be argued that the first ACF minimum occurs at lag = and the first ACF maximum 
occurs at lag = P^. On the other hand, when ^ successive minima and maxima in the ACF 
are separated by the period P^ of the spike train. In this regime, the ACF carries no information 
about the width of a spike, and we expect that the lag at the first maximum (be., P^) = 2 x the 
lag at the first minimum (be.,t^). Conversely, any indication of the relationship P^ = 2 x in 
the inferred timescales may therefore be interpreted as suggesting the lack of complete information 
about the two timescales. 


To motivate the methodology we are about to develop for inferring microstructnre timescales 
from individual subpnlses in a pulsar data set, consider four illustrative subpnlses from PSR 
B0525+21 together with their ACFs (Fig. 4). In the left column of this figure, black curves repre¬ 
sent probable candidates for the signal in the data time series (gray scatter), red curves represent 
probable candidates for the smooth trend (be., envelope) in the pulse, and bine curves represent 
probable candidates for the microstructnre (obtained after removing the probable envelope from 
the probable signal). In the right column of this figure, the gray curves represent the ACF for the 
whole-pulse time series, whereas the bine curves represent the ACF for the microstructnre feature. 
Note that while the ACF for the whole-pulse time series must be non-negative at all lags (because 
pulse time series are by-and-large positive, barring noise), the ACF for the microstructure feature 
may be negative (because the microstrncture signal fiuctuates on both sides of 0). In Fig. 4, the 
time series for pulse 179 shows strong microstructnre oscillations superimposed on top of a smooth 
envelope, which translate into weak but clear oscillations in the descending part of the snbpulse 
ACF close to lag 0. For this pulse, we see that the first minimum in the subpulse ACF (gray curve 
in the ACF plot) agrees with the first minimum in the ACF for the microstructnre feature (blue 
curve in the ACF plot). This pulse is an example where the conventional method of identifying 
microstrncture timescales using the first minimum in the ACF will work flawlessly because the first 
minimum in the subpulse ACF is numerically well-defined. The time series for pulse 37 again shows 
a highly finctnating signal characteristic of the microstrncture. However, the first ACF minimum 
is barely visible in the snbpulse ACF. In contrast, the first minimum is clearly identifiable in the 
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microstructure ACF for this pulse. The time series for pulse 21 again shows a clearly identifiable 
microstructure feature. However, it nevertheless leads to a featureless subpulse ACF which makes 
it impossible for the conventional method to identify a microstructure timescale from it. The ACF 
for the microstructure feature, in contrast, has a clearly identifiable first minimum. Pulse 62 is an 
example where, with the exception of one or two clear but weak fiuctnations in the time series, 
there may not be any significant/detectable microstrnctnre feature. 

These illustrative pulses suggest that the signal hidden in the data can be modeled as con¬ 
sisting of an envelope feature characterized by relatively strong long-time smooth variations, and 
a micro structure feature characterized by relatively weak short-time qnasiperiodic variations that 
ride on top of the envelope. Hence, with reference to the use of ACFs for identifying timescales, 
we expect the following two confounding situations: 

1 . The ACF of the complete pulse may be dominated by power at low frequencies, be., by 
the envelope feature, possibly masking any short-time modulations (e.^., pulses 37 and 36 in 
Fig. 4). 

2 . Noise in the data may result in additional random oscillations in the ACF, making it nu¬ 
merically difficult to identify genuine minima, slope changes, etc., in the ACF for the whole 
subpulse. 

This suggests that if the noise and the envelope feature are extracted from the data to arrive at 
the microstrnctnre feature, then the microstrnctnre feature can be used meaningfully to estimate 
microstrnctnre timescales from its ACF in the conventional manner. 

Therefore, a solution to the problem of estimating microstrnctnre timescales lies in correctly 
estimating and removing the noise as well as the long-time smooth modulation (be., envelope) from 
a snbpulse time series. Our methodology for estimating microstrnctnre timescales from individual 
subpnlses thus consists of the following steps. 

1. The average pulse profiles for many pulsars are often seen to be composed of multiple emit¬ 
ting components. As was pointed out in some earlier studies (e.^-., Soglasnov et al. 1981; 
Soglasnov et al. 1983; Popov et al. 2002), the individual components can have different mi¬ 
crostructure timescales. Hence, we isolate individual pulse components wherever possible 
before estimating timescales. Individual components are identified from the average pulse 
profile in the form of time windows (shown as gray sections in the figures in Sec. B) relative 
to the beginning of the profile. These windows are then applied to each single pulse in the 
data set to separate it into its components. For pulsars where multiple components could not 
be identified, we used whole pulses. Further, we restrict our analysis only to pulses with high 
S/N, where S/N is estimated as the ratio of the peak signal to the off-pulse noise. In Ta¬ 
ble 1, 24 pulsars marked in bold face are the ones with sufficient single-pulse S/N and hence 
were suitable for timescale analysis, while the others are not used. The pulsar rotational 
parameters for these 24 pulsars are given in Table 2. 
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2 . RFI-induced outlier values in subpulse time series data may influence timescale estimates. 
We therefore identify potential outliers in subpulse time series data using a heuristic method, 
and replace them with appropriate values resembling their respective local neighborhoods in 
the time series (Sec. A.2). The decision to “curate” outliers in this fashion is made on a 
case-to-case basis through preliminary visual inspection of each pulsar data set. 

3. We de-noise every snbpulse using a model-independent smoothing spline fit (Sec. A.3) that 
adapts to the shape of the subpnlse optimally. The fit thus obtained is taken to be a smooth 
noise-free version of the snbpulse data (e.^-.^the black curves in Fig. 4 left column). 

4. We estimate the envelope for each snbpulse via kernel regression (Sec. A.4) using a heuristic 
smoothing timescale. The difference between the fit and this envelope is the estimated mi- 
crostrncture. Note that the smoothing timescale (referred to as smoothing bandwidth h later 
on) is the most critical parameter for our methodology. Further details can be found in Sec. 

4. 

5. The nature and strength of the microstrnctnre feature can vary considerably among snbpulses. 
Specifically, the microstructure feature may be nearly absent or at best weak in some of the 
snbpulses. We identify subpulses with weak microstructure using the degrees of freedom (DF) 
for the smoothing spline fit (Sec. A.3), because lower DF values correspond to smoother fits 
and weaker microstructures. For better estimation of microstructure timescales, we reject 
subpulses with DF in the lower 10-20 percentiles of the DF distribution for that data set (See 
Supplementary Table SI for the percentile cut-off for each data set). 

6 . For the selected pulses, the two characteristic microstrnctnre timescales (width and pe¬ 
riod Pfj) are estimated respectively as the first minimum and the first maximum in the ACF 
of the estimated microstructure for a snbpulse. Because of the substantial variation in mi¬ 
crostructure characteristics across snbpulses, the timescales thus obtained need to be further 
summarized or analyzed statistically. Specifically, for each pulsar, we report the median 
timescale together with the median absolute deviation (MAD) as a measure of spread. We 
also explore the nature of any relationship between and P^. 

The microstructure timescale analysis pipeline developed in this paper is implemented in the R 
(R Core Team 2014) open source statistical computing environment. Detailed pnlsarwise results is 
available as a tarball in: ftp://wm.ncra.tifr.res.in/dmitra/hires_analysis_pdf.tar.gz or 
http://cms.unipune.ac.in/reports/tr-20150223/. 


4. Results 

Fig. 5 and 6 depict typical results from our analysis package. Fig. 5 shows how the technique 
is applied to a snbpulse (see caption for description). We estimate the microstrnctnre periodicity 



as the first peak in the ACF of the microstructure signal (blue line), and show the results 
for three envelope smoothing bandwidths h because the estimated timescales are sensitive to the 
envelope bandwidth, where h is the smoothing window size in units of the total number of points in 
a subpnlse (see Sec. A.4 for details). In the absence of any theoretical prediction for microstrnctnre 
timescales, it is difficult to choose an appropriate value for h from first principles. Hence, we used 
timescale estimates from earlier microstructure studies in the literature to guide us, together with 
crucial properties of the subpulse data to converge upon a range of reasonable h values, as follows. 
For every pulsar, we identified a few highly fiuctnating/modulating pulses (like examples shown 
in Fig. 1) where the P^ estimated from first peak in the ACF is easily discernible. This gave ns 
a range of h values over which the P^ value remained unchanged. Using a few representative h 
values in this range we found the distributions of P^. The median value of these P^ distribution 
(or found from first dip in the ACF) appeared to be in good agreement with available estimates 
in the literature. We found that for all pulsars, h larger than 0.1 produced envelopes that were too 
smooth, whereas h less than 0.05 led to envelopes that traced the subpnlse data too closely. All 
our results are hence quoted for three representative bandwidths, namely, 0.1, 0.075 and 0.05, as 
illustrated in Fig. 6. The median of the P^ histogram is given in columns 5, 6 and 7 of Table 3 
where the spread of the P^ distribution is estimated using the median absolute deviation (MAD). 
Column 2 gives the longitude range used to define the subpnlses for microstrnctnre analysis, and 
these can be identified from the average pulse profiles presented in Sec. B. 

The top panel of each plot in Fig. 6 deserves attention. The plot depicts the relationship 
between the values of the first dip (t^) along the z-axis and the first peak (P^) along the y- 
axis as estimated from the ACFs of the microstrnctnre. The three dotted lines correspond to 
= 2tf^ and P^ = respectively. Most of the points in this figure (and also for 
other pulsars) are seen to cluster around the P^ = 2t^ line, indicating that the microstrnctnre is 
unresolved and hence the width cannot be determined. 

Before we close this section, it is important to compare our analysis method with earlier 
studies. Majority of the studies use the first dip or peak of the ACF to estimate microstrnctnre 
timescales (e.^., Soglasnov et al. 1981; Hankins 1972). However, as mentioned earlier, the broad 
subpnlse power as well as noise in the subpnlse make it difficult to extract the timescales nu¬ 
merically. Our method is designed to alleviate these problems by first de-noising the subpnlse 
profile and then removing a smooth envelope from it. Other studies have used the power spec¬ 
trum to estimate timescales corresponding to quasiperiodic microstructures {e.g., Popov et al. 2002; 
Lange et al. 1998). The slow-varying envelope in these cases appears as low-frequency power in the 
spectrum. For example, Popov et al. (2002) constructs the ACF by cross-correlating the voltages 
from two different spectral channels, while Lange et al. (1998) use the derivative of the power spec¬ 
trum, although noise in single pulses may still affect their analyses. The strength of our method 
lies in its ability to recover low-amplitude intensity modulations in a subpnlse, isolating the mi¬ 
crostrnctnre from the subpnlse time series. Consequently, we find that majority of the subpnlses 
have microstructure emission superposed on a smooth broad emission (envelope). 
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5. Distribution of timescales for all Stokes parameters and the P^-P relation 

Cordes (1979), based on 10 pulsars observed at 430 MHz in the pulsar period range of 0.2 to 3.7 
second, found the relationship t^ ~ 10 “^P between the microstrnctnre timescale and the pulsar 
period P. Kramer et al. ( 2002 ) reestablished this relationship as t^(/is) ~ (600 ± 100)P(s)^-^^°-^ 
using 12 pulsars where estimates were available in the literature for different frequencies. With 
our current set of observations we are in a position to revisit this relationship using 24 pulsars at L 
Band. Since broad microstructures cannot be resolved, we believe that the periodicity P^ is a more 
relevant quantity instead of Fig. 7 plots the the median value of P^ against period P. The top 
panel corresponds to the total intensity, the middle one to the linear polarization, and the bottom 
one to the circular. The P^ points are estimated from the median value of the histogram of the first 
peak in the ACF, and the error bars reflect the median absolute deviation of the P^ distribution. 
A linear fit of the form P^ (msec) = mP (msec) + c was performed using the median P^ values 
for stokes I corresponding to h = 0.1,0.075 and 0.05 separately, which yielded (m, c) values as 
(1.3 X 10“^, 0.04), (1.1 X 10“^, 0.08) and (0.9 x 10“^, 0.098) respectively. These fits are shown in 
Fig. 7. We note that the fit is merely demonstrating the tendency that median P^ increases with 
increasing pulsar period P, although a large spread of P^’s exists for all pulsars. 

Correlation of P^ with other pulsar rotational parameters like period derivative, characteristic 
age and surface magnetic field did not show any significant trends. Although a weak anticorrelation 
was seen (not shown here) between the slowdown energy (P) and P^, more data is needed to 
establish this. 


6. How does broad microstructure average to form a subpulse? 

It was evident from early polarization observations (e.g., Clark & Smith 1969; Manchester et al. 
1975) that single pulses are highly linearly polarized, and that the circular polarization changes sign 
or handedness within a pulse. Based on high-quality single-pulse polarimetry, Mitra et al. (2009, 
hereafter MGM09) have showcased a collection of snbpulses from an assorted set of pulsars, where 
the subpulses had a Gaussian shape, with close to 100% linear polarization and associated with 
sign-changing circular. One such example is shown in the left panel of Fig. 8 . 

For a long time this emission pattern kindled great interest in the pulsar community (see 
e. 5 -., Michel 1987; Radhakrishnan & Rankin 1990), because this remarkably resembles the emission 
pattern due to single-particle vacuum curvature radiation as the observer cuts across the emission 
cone {i.e., the I /7 cone). The mathematical formulation for obtaining Stokes parameters due to cur¬ 
vature radiation in vacuum for a single relativistically moving particle has been explicitly calculated 
by several authors (e.^.,Gil & Snakowski 1990b; Gil k, Snakowski 1990a; Ahmadi k Gangadhara 
2002 ). Let ns recall the polarization properties of single-particle curvature radiation in vacuum 
based on the description given by Gil k Snakowski (1990a), Gil k Snakowski (1990b). Gonsider 
a charged particle moving relativistically with Lorentz factor 7 , along super-strong magnetic field 
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lines, with a radius of curvature p. The particle will emit curvature radiation beamed in the forward 
direction. The Fourier components of the electric fields can be decomposed into parallel (£||) and 
perpendicular components to the plane of the curved magnetic field lines (see eqs. 7 and 8 of 
Gil & Snakowski 1990a). When the radiation is viewed as a cut through the emission beam, the 
has a sinusoidal form, £|| is positive and has a relative phase-shift of 90° forming an orthogonal 
basis. The radiation polarized in the parallel and perpendicular direction are also known as the 
ordinary (O) and extraordinary (X) modes. The power in curvature radiation is dependent on the 
frequency of emission, and achieves its maximum at tfie characteristic frequency which is given 
by r'c = l.SSy^r'o fbe vacuum case (See fig. 2 of Gil & Snakowski 1990b). The ratio of the two 
components as well as the width of the pulse is dependent both on the viewing angle (f and the 
frequency of emission (see figs. 3, 4 and 5 of Gil & Snakowski 1990b). The width of the pulse is 
approximately equal to 4 = I /7 at Uc- The power in the O mode is seven times stronger than the 
X mode in vacuum (Jackson 1975; Gil et al. 2004). 

Our aim here is to look for signatures of the aforementioned vacuum curvature radiation in 
the data. First, let ns note that the subpulse width cannot be associated with the elementary 
1/7 emission cone, because when observed with higher time resolution, the subpulse splits into 
microstructnres. In an attempt to reproduce this snbpnlse radiation pattern, Gil & Snakowski 
(1990a) and Gil et al. (1993) used the shot-noise model of pulsar emission, where the subpulse is 
an incoherent addition of smaller bunches, with each bunch emitting like a single-particle vacuum 
curvature radiation. The subpulse shape results from the intensity gradient due to the charge den¬ 
sity variation of the bunches across the subpulse. These authors demonstrated that in the presence 
of intensity gradients, the resulting emission will have a net circular polarization as observed in 
snbpulses. The beams from the bunches do not overlap, so the linear polarization remains high. 
Assuming 7 ~ 100, and therefore an angular size 1 / 7 ^ 0.6° for the emitting bunches, they simu¬ 
lated several pulses that appear similar to the snbpulses showcased by MGM09. These results were 
also recently confirmed by Gangadhara ( 2010 ). A prediction of this model is that if a subpulse 
can be resolved into microstructnres (e.^-.^see fig. 3 of Gil et al. 1993), then the microstructnre 
should have features which mimic the single-particle vacuum curvature-radiation pattern. To test 
this theory we need examples of snbpulses where the microstructnres are highly resolved. Note 
that the snbpulses chosen should be highly linearly polarized such that they do not suffer from any 
depolarization effects due to orthogonal polarization moding^ or due to incoherent addition from 
overlapping I /7 emission beam (see MGM09). 

In the left hand panel of Fig. 8 , a subpulse of PSR B2020+28 is shown which is highly linearly 


^It is important to touch upon a subtle point regarding sign-changing circular polarization in pulsar subpulses. 
In some cases, at the zero-crossing point of the sign-changing circular, the linear polarization drops significantly and 
the polarization position angle (PPA) gets associated with an orthogonal jump. However, for testing the vacuum 
curvature-radiation signatures, such subpulses are undesirable as they do not comply with the basic property of 
vacuum curvature radiation where, at the zero-crossing point of the circular, the fractional linear polarization has to 
be close to 100%. 
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polarized and is associated with a sign-changing circular polarization, similar to the pulses show¬ 
cased by MGM09. The right hand panel shows the high time resolution data for the same subpulse, 
and clearly resolved microstructures are seen. As expected, the linear polarization of the individual 
microstructures are close to 100% in order to produce the highly polarized snbpulse. The circular 
polarization below the microstructnre feature does not show any sign reversals. Rather it has the 
same width (and also periodicity) as the total intensity and linear polarization. In fact, this de¬ 
scription is true for all snbpulses where sufficiently resolved microstructnres are seen. Fig. 1 shows 
several such examples, and very clear examples where Stokes V has the same as Stokes I and L 
are shown in Fig. 9. This feature is also borne out from the statistical analysis of the data set where 
the median values of Stokes /, L and V are similar as can be seen in Table 3. Thus, based on 
these observations, one is compelled to conclude that radiation pattern of resolved microstructnres 
are not consistent with single-particle curvature radiation of charged bunches in vacuum. In the 
next section, we will discuss how these findings corroborate our current understanding of the pulsar 
emission mechanism. 


7. Pulsar microstructure and its implications for the coherent radio emission 

mechanism 

Coherent pulsar radio emission can be generated by means of either a maser or coherent 
curvature mechanism (Ginzburg & Zhelezniakov 1975; Kazbegi et al. 1991; Melikidze et al. 2000, 
hereafter MGPOO). There is general agreement that this radiation is emitted in curved magnetic 
field-line planes due to growth of plasma instabilities in the strongly magnetized electron-positron 
plasma well inside the light cylinder. Many observational constraints on the emission altitude 
suggest that the emitted radiation detaches from the ambient plasma at altitudes of about 500 km 
above the stellar surface which is typically less than 10% of the light cylinder radius P£,c = Pcj2 ti 
(e.^.,Cordes 1978; Blaskiewicz et al. 1991; Rankin 1993; Kijak & Gil 1997; Mitra & Deshpande 
1999; Gangadhara & Gupta 2001; Mitra h Li 2004; Krzeszowski et al. 2009). The magnetic field 
at these altitudes is so strong that the plasma is constrained to move only along the magnetic 
field lines. Under such conditions, the two-stream instability is the only known instability that 
can develop, and the growth of this plasma instability leads to the formation of charged relativistic 
solitons, which can excite coherent curvature radiation in the plasma {e.g., Asseo & Melikidze 1998; 
Melikidze et al. 2000; Gil et al. 2004 GLM04 hereafter). 

The basic requirement of the theory of coherent radio emission is the presence of a non¬ 
stationary inner acceleration region. Ruderman & Sutherland (1975) (hereafter RS75) were amongst 
the first to explore the conditions in the pulsar magnetosphere that lead to the formation of an 
inner accelerating region (also known as an inner vacuum gap, IVG) just above the pulsar polar 
cap where the free fiow of charge is hindered, thus resulting in the formation of an inner vacuum 
gap with a large potential difierence. The physics of inner vacuum gap formation can be found in 
RS75, and here we point out that the original RS75 inner vacuum gap was not consistent with x-ray 
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observations and radio subpulse drifting observations, and currently it is thought that a partially 
screened inner vacuum gap exists close to the neutron star surface (Gil et al. 2003). The process 
that generates non-stationary flow of charges in the gap is as follows. The IVG is initially dis¬ 
charged by a photon-induced electron-positron pair-creation process in the strong curved magnetic 
held. The height (hg) of the gap stabilizes typically at about a few mean-free paths of the pair- 
production process. The electric held in the gap causes separation of the electrons and positrons 
(also known as primary particles), where the electrons bombard onto the neutron-star surface, and 
the positrons are accelerated away toward the outer magnetosphere. As these primary particles 
move they accelerate and create further secondary particles resulting in a pair-creation cascade. 
The outflowing particles are hence composed of primary particles with Lorentz factors of about 
10® and secondary plasma clouds with a distribution of Lorentz factors from about 10 to 1000, but 
peaking around 300-500. This process eventually occurs over a set of magnetic held lines, and this 
overall spread of the discharging process is called a spark. As the electron moves towards the polar 
cap it experiences E x B drift, causing the sparking discharge to move slowly across the polar cap. 
RS75 envisaged that the inner vacuum gap is discharged in the form closely packed isolated sparks, 
where the lateral size and distance between the sparks is equal to the gap height. A plasma column 
is associated with each spark, and in the radio-emission region coherent spark-associated plasma 
columns radiate and generate the observed subpulses in pulsars. The sparking process continues 
for a time G until the whole potential in the gap is screened out, after which the plasma empties 
the gap in a time tg = hg/c and once again the sparking process starts almost in the same place 
on a timescale ^ tg « tg (Gil & Sendyk 2000). This process results in a non-stationary flow of 
plasma clouds associated with every spark anchored to a speciflc location on the polar cap. The 
physics of the sparking discharge process is still a topic of active research, however most studies 
suggest that in normal pulsars hg is of the order of about 50-100 m, and to screen the gap potential 
around 50 pair creations are needed; and hence G ^ hOhg/c ~ 10// sec. Note that during the time 
G the E X B drift causes the spark to move only very slightly across the polar cap. Two successive 
plasma clouds, with each cloud extending for about hOhg and separated by a distance hg, even¬ 
tually overlap when the fast moving particles of the second cloud catch up with the slow moving 
particles of the first cloud, leading to the development of the two-stream instability. The instability 
in the plasma generates Langmuir plasma waves and the modulational instability of the Langmuir 
waves leads to the formation of a large number of charged solitons which can excite extraordinary 
(X) and ordinary (O) modes of curvature radiation in the plasma cloud as shown by MGPOO and 
GLM04. The condition for coherent emission requires the wavelength of emission to be greater than 
the intrinsic soliton size, < 2y/^ujp where ujp = (dvre^n/me)^ is the plasma frequency (see 

MGPOO, GLM04), and is the Lorentz factor of the secondary plasma. The emerging radiation 
is an incoherent addition of the radiation emitted by each soliton. 

We now discuss aspects of the soliton coherent-curvature radiation theory that relate to the two 
important microstrncture behaviors discussed in this paper, namely the periodic structures observed 
in pulsar subpulses, their linear dependence with pulsar period and the absence of signatures of 
vacuum curvature radiation in pulsar microstrnctures. Obviously the timescales involved in the 
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radiation process are far too small to be seen in our 59.5 //sec observations, hence we rule out that 
the periodic microstrnctures correspond to emission from individual plasma clouds or solitons^. 
Another way of sampling a periodic structure would be if temporal modulation of the emission 
exists in the emission region. There are theoretical studies in the literature which attempt to 
interpret microstructure based on nonlinear plasma effects leading to temporal modulation of the 
pulsar emission (e.^., Chian & Kennel 1983; Weatherall 1998; Asseo 1993); however these models 
are mostly developed for emission where uj > 2^/^ujp, whereas coherent emission requires ou < oup. 

The other important aspect of our observations is the absence of a signature of vacuum cur¬ 
vature radiation in pulsar microstrnctures. The coherent curvature radiation by charged bunches 
called solitons excites the extraordinary and ordinary mode in the electron-positron plasma. Due 
to the strong ambient magnetic field the X-mode does not interact with the plasma and hence 
can escape the plasma and reach the observer (MGM09, Melikidze et al. 2014) . The fate of the 
0-mode is however unclear. It has been suggested that it can be damped and gets ducted along the 
magnetic field lines, or if there are sharp plasma boundaries that it can escape as electromagnetic 
waves and reach the observer. In any case, due to propagation effects in the plasma, the X and O 
modes separate in phase and hence they do not maintain any phase relation as they emerge from 
the plasma. As a result the circular polarization can no longer have any sign-changing signature 
as in the vacuum case. Also incoherent addition of a large number of solitons would tend to cancel 
any circular polarization, particularly in our observations. Thus the origin of circular polarization 
still remains to be understood in this model of pulsar emission. 


8. Summary 

In this work we have found several distinguishing features of pulsar microstrnctures. These 
can be summarized as follows 

1 . Almost all subpulses show a qnasiperiodic structure riding on top of a broad pulse envelope. 
In general this structure is weak, except in occasional cases when deep modulations are seen. 
However, the absence of microstrncture in pulsar average profiles is a clear indication that at 
any given longitude pulsar emission is random and stochastic. 

2 . Microstrnctures are also seen to be associated with both the orthogonal polarization modes. 

3. We did not find any case where microstrncture widths are resolved. In other words, we find 
that Pp ^ 2tp (where the estimated width is tp found as the first dip in the ACT) is related 
to the periodicity Pp (found as the first peak in the ACT) (see Sec. 3 for details). 


^Note that in this context the findings of Popov et al. (2002) are rather interesting as based on their 62.5 nanosec 
observations they report that the shortest microstrnctnre widths are abont 2/xs to 10/xs, which corresponds to the 
time over which the sparking discharge takes place. 
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4. For every pulsar the values extracted from a collection of subpulses have a distribution, 
and the median value of linearly increases with pulsar period (see Sec. 5). 

5. All the Stokes parameters have the same periodicity—be.,they have same values (see 
Sec 6). 

A variety of analysis techniques to extract microstructure timescales have been deployed in 
the past. The methods use either the ACF first dip and first peak or use the peaks in the power 
spectra for estimating timescales. However, these methods cannot effectively find the low-amplitude 
wiggles that are observed in snbpulses. We address this problem by de-noising the subpulse signal 
and subsequently removing a broad envelope. The ACF of the resultant signal is then used to 
estimate microstructure timescales. 

These qnasiperiodic structures, which show a linear dependence on rotation period, cannot be 
understood in terms of the current theories of pulsar emission. We also find that the linear and 
circular polarization follow the same distribution of timescales as the total intensity. This, as we 
have argued strongly, rules out the possibility that microstrnctnres are basic units of curvature 
radiation in vacuum. We conclude that the microstrnctnres within the pulsar signal reflect tem¬ 
poral modulations of plasma processes of the emission engine. While this paper is an exercise in 
understanding the properties of pulsar microstrnctnres, theoretical progress needs to be made to 
understand this unique phenomenon. 
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Fig. 1.— A collection of subpulses are shown in full polarization for various of the pulsars observed 
with a time resolution of 59.5 fisec. In each plot the top panel shows Stokes /, L and V in black, 
red and blue lines respectively. The polarization position angle (PPA) histogram of the entire pulse 
sequence is shown as a “dotty plot” display of qualifying samples in the bottom panel. The average 
PPA of the subpulse is plotted in green when the circular polarization is negative and magenta 
when positive. The ringing/modnlating structures seen in the subpulses are microstructures. Note 
that a certain periodicity (or qnasiperiodicity) is apparent in the microstructures, and all the Stokes 
parameters fluctuate in the same manner with the same width and periodicity. 
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B1133+16 Aracibo 1177.66 MHz 1.186 b IIIPuIm B1133-I-16 Arecibo 1177.66 MHz 1.168 a 135Pljl8e 
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Fig. 2.— Two single pulses of PSR B1133+16 (plot description same as given in the caption of 
Fig. 1), illustrating that microstructures occur within both polarization modes. In the left hand 
display around 0° longitude the pulse component showing the microstructure feature corresponds 
to the weaker secondary polarization mode; whereas in the right hand one at about 6° the mi- 
crostructnre corresponds to the primary polarization mode. 
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Fig. 3.— A toy example illustrating the relationship between the width of a spike, period of a 
spike train, and the first minimum and maximum in the ACF. In the top row, each spike is a 
truncated Gaussian with baseline width = 200, and the spike train is a sequence of equidistant 
spikes separated with the same gap. The period of the spike train therefore equals the widths of 
a spike + gap between successive spikes, (a) width < gap; (b) width = gap; (c) width > gap. In 
all cases, the first ACF maximum is at lag = period. For (a) and (b), the first ACF minimum is at 
lag = width, and the ACF can be used to infer both the width and period of the spike train. For 
(c), the first ACF minimum is at lag = period/2, and the ACF carries no information about the 
width. The spike train is a caricature of the microstructure feature of a subpulse, and illustrates 
what can be inferred about timescales using the ACF. Note that the bottom row plots logio(ACF) 
instead of the ACF so as to accentuate the ACF minima. 
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PSR B0525+21 | Pulse 179 



0.71 ms 
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PSR B0525+21 | Pulse 37 




PSR B0525+21 | Pulse 36 




PSR B0525+21 | Pulse 62 




Fig. 4.— A few illustrative subpulses from a B0525+21 data set and their respective ACFs. Left 
column: gray, subpulse time series data; black, probable signal in the data sans the noise; red, 
probable smooth trend in the time series; blue, probable microstructure feature of the signal. Right 
column: gray, ACF computed directly from noisy data, blue: ACF computed from the probable 
microstructure feature. See text for detailed description. 
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Fig. 5.— The analysis method for extracting the timescale of a single subpulse for the trailing 
component of PSR B1133+16 is illustrated in this figure. The plot has three columns, with each 
column corresponding to Stokes I, L and V, respectively. In each column the upper panel shows 
the snbpulse, where the gray line corresponds to the snbpulse data, the black line is the fit to the 
data, the red line is the envelope for a given bandwidth and the blue line is the microstructnre 
signal defined as the difference between the red and black lines. The second panel correspond to 
the ACF of the microstructnre (blue) signal. The timescales that we measure are indicated in the 
figure. The three sets of plots from top to bottom are for three different bandwidth values 0.1, 
0.075 and 0.05, respectively (see text for details). 
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Envelope smoothing bandwidth: 0.075 
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Fig. 6.— The three plots above are histograms of the estimated timescales of the microstructure 
signal for the trailing component of PSR B1133+16. The three sets of plots correspond to three 
different bandwidths h = 0.1,0.075 and 0.05, respectively, from top to bottom . In each plot there 
are three columns corresponding to Stokes I, L and V from left to right, respectively. The topmost 
panel in each plot is the relation between the first dip (as x-axis) and first peak (as y-axis) of the 
ACF. The middle panel shows distribution of the first dip in the ACF which corresponds to 
The bottom panel shows the distribution of the first peak in the ACF which correspond to (see 
text for details). 
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versus Pulsar period (Stokes I) 



Pulsar Period, P [msec] 

P^ versus Pulsar period (Stokes L) 



Pulsar Period, P [msec] 

P^ versus Pulsar period (Stokes V) 



Fig. 7.— The filled points and the error bars in the above plots are the median value and the 
median absolute deviation of distribution given in Table 3 for all the pulsars as a function of 
pulsar period P. The three plots from top to bottom correspond to Stokes I, L and V, respectively. 
The red, green and blue points corresponds to bandwidths of 0.1, 0.075 and 0.05, respectively. 
The dashed red, green and blue line is a fit to the median P^ values for stokes I corresponding 
to bandwidth 0.1, 0.075 and 0.05, respectively. The same lines are shown for stoke L and V for 
reference. The orange line at 300 jusec along the abscissa is about five times the time resolution of 
our observations. 
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B2020+28 Arecibo 1399.92 MHz 0.343 s 125Pulse 


B2020+28 Arecibo 1399.92 MHz 0.343 s 125Pulse 



Longitude [ ° ] 


Longitude [ ° ] 


Fig. 8.— The left plot shows a subpulse of PSR B2020+28 (plot description same as given in the 
caption of Fig. 1), which has been smoothed to time resolution of 595 yusec. Note that the subpulse 
is close to 100% linearly polarized and shows a sign-changing circular. The corresponding right 
display shows the high resolution pulse with time resolution 59.5 /rsec. The striking point to note 
there is that the subpulse is split into microstructures, with the linear and circular polarization 
following the same width and periodicity as that of the total intensity. These microstructures hence 
do not represent coherent curvature radiation by charged bunches in vacuum (see text for details). 
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Fig. 9.— A few archetypal examples (as in Fig. 1) illustrating that the circular polarization follows 
the same values as the total intensity and linear polarization. In each display a highly resolved 
subpulse is shown, and in the bottom panel the ACF of the Stokes parameters are shown. Note that 
in all the examples the PPA across the snbpulse is primarily associated with a single polarization 
mode. Clearly the ACFs of Stokes I, L and V trace each other extremely well. 
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Table 1. Observational Parameters: Column 1 to 3 in the table below lists the pulsar name, 
observing frequency and the observing MJD. Column 4 gives the channel bandwidth and the total 
number of channels across the band. Column 5 gives the total number of pulses 
observed/resolution in degrees/resolntion in //sec. Column 6 and 7 gives the pulsar rotation 
measure and the dispersion measure (DM) used for our analysis. Column 8 and 9 gives the pulse 
averaged percentage linear and circular polarization respectively. Note that pulsar name with 
boldface are the ones used for microstrnctnre timescale estimate analysis. 


PSR 

Name 

Band 

(MHz) 

MJD 

Ch.BW (GHz)/ 
Nchan 

length 
res (°)//xsec 

RM 

rad 

DM 

pc/cc 

%L 

%V 

B0301+19 

337.668 

55283 

0.085/256 

616/0.015/59.52 

-8.3 

15.74 

43 

17 

B0525+21 

351.764 

56572 

0.085/256 

616/0.015/59.52 

-39.6 

50.94 

24 

-1.9 

J0546+2441 

351.764 

56573 

0.024/1024 

843/0.021/166.67 

-19.6 

73.81 

6 

-3.2 

B1919+21 

352.806 

55320 

0.112/512 

1139/0.016/62.5 

-16.5 

12.46 

14 

8 

B1929+10 

337.668 

55320 

0.085/256 

4639/0.094/59.52 

-5.6 

3.18 

88 

-3.6 

B1944+17 

352.806 

55320 

0.112/512 

2269/0.051/62.5 

-28.0 

16.22 

6 

-6.6 

B2016+28 

337.668 

55320 

0.085/256 

953/0.038/59.52 

-34.6 

14.17 

23 

-0.7 

B2020+28 

337.668 

55320 

0.085/256 

4211/0.062/59.52 

-74.7 

24.64 

44 

0.5 

B2044+15 

337.668 

55320 

0.085/256 

960/0.019/59.52 

-100.0 

39.84 

36 

-10.3 

B2110+27 

337.668 

55320 

0.085/256 

1265/0.018/59.52 

-37.0 

25.11 

4 

-3.3 

B2315+21 

352.806 

55275 

0.112/512 

275/0.015/59.52 

-35.8 

20.91 

32 

3 

B0301+19 

1220.668 

55425 

0.085/1024 

1202/0.015/59.52 

-8.3 

15.74 

36 

8.8 

B0525+21 

1220.668 

55522 

0.085/1024 

640/0.006/59.52 

-39.6 

50.94 

45 

-13.7 

J0546+2441 

1220.544 

56572 

0.085/1024 

421/0.007/59.52 

-19.2 

73.81 

35 

-13 

B0656+14 

1220.668 

55522 

0.085/1024 

1558/0.055/59.52 

-23.5 

13.997 

91 

-19.7 

B0751+32 

1220.668 

55522 

0.085/1024 

1944/0.014/59.52 

-7.0 

39.949 

29 

-26 

B0823+26 

1220.668 

55522 

0.085/1024 

1130/0.040/59.52 

5.9 

19.454 

19 

0 

B0834+06 

1220.668 

55522 

0.085/1024 

1006/0.017/59.52 

23.6 

12.889 

15 

-4.5 

B0919+06 

1220.668 

55522 

0.085/1024 

1045/0.049/59.52 

29.2 

27.271 

48 

4.6 

B0950+08 

1220.668 

55522 

0.085/1024 

2055/0.084/59.52 

0.66 

2.958 

22 

-4.5 

B1133+16 

1220.668 

55522 

0.085/1024 

2104/0.018/59.52 

1.1 

4.864 

28 

-9 

B1237+25 

1220.668 

56550 

0.085/1024 

873/0.015/59.52 

-0.33 

9.242 

51 

9 

J1503+2111 

1220.416 

56550 

0.085/1024 

1025/0.006/59.52 

16.0 

11.75 

22 

30 

B1534+12 

1663.008 

55637 

0.067/512 

15835/0.36/38.69 

10.6 

11.614 

30 

-2.9 

J1713+0747 

1220.668 

55632 

0.085/1024 

65646/4.735/59.52 

8.4 

15.9915 

35 

-1.4 

J1720+2150 

1220.544 

56550 

0.086/1024 

1021/0.013/59.52 

48.0±10 

41.4 

39 

-63 

B1737+13 

1220.668 

55632 

0.085/1024 

747/0.026/59.52 

64.4 

48.673 

45 

-3.2 

J1740+1000 

1220.668 

55632 

0.085/1024 

3893/0.139/59.52 

23.8 

23.85 

88 

-23 

J1746+2245 

1220.544 

56550 

0.085/1024 

899/0.006/59.52 

91.0±5 

52.0 

29 

42 

B1855+09 

1220.668 

55637 

0.086/1024 

111910/4./59.52 

16.4 

13.295 

22 

1.8 

B1901+10 

1220.752 

56563 

0.085/1024 

1034/0.011/59.52 

-105±5 

135 

36 

35 

J1910+0714 

1220.752 

56563 

0.085/1024 

1024/0.007/59.52 

182±20 

124.06 

20 

-12 

B1910+20 

1220.752 

56563 

0.085/1024 

1025/0.009/59.52 

148±10 

88.34 

37 

-17 

B1919+21 

1563.008 

56207 

0.672/128 

346/0.016/59.52 

-16.5 

12.46 

10 

0.6 

B1929+10 

1220.668 

55633 

0.085/1024 

2648/0.095/59.52 

5.26 

3.180 

68 

-26 

B1944+17 

1220.668 

55633 

0.085/1024 

1361/0.048/59.52 

28.0 

16.22 

22 

-14 

B2002+31 

1220.668 

56564 

0.085/1024 

1054/0.010/59.52 

30.0 

234.82 

15 

2.6 



Table 1 


Continued 


PSR 

Name 

Band 

(MHz) 

MJD 

Ch.BW (GHz)/ 
Nchan 

length 
res (°)//isec 

RM 

rad 

DM 

pc/cc 

%L 

%V 

B2016+28 

1220.668 

55632 

0.085/1024 

1075/0.038/59.52 

-34.6 

14.172 

23 

-4 

B2020+28 

1220.668 

55632 

0.085/1024 

1747/0.062/59.52 

-74.7 

24.64 

37 

-6.5 

B2034+19 

1220.627 

56564 

0.085/1024 

1341/0.014/59.52 

-97±10 

36.0 

24 

5 

B2110+27 

1470.668 

55425 

0.085/1024 

1625/0.018/59.52 

-37.0 

25.113 

13 

15 

B2113+14 

1470.668 

55425 

0.085/1024 

1060/0.048/59.52 

-25.0 

56.149 

19 

-15 

B2315+21 

1470.668 

55425 

0.085/1024 

746/0.014/59.52 

-35.0 

20.906 

17 

-15 



Table 2. Pulsar Parameters used for microstructure study in this paper obtained from the 
ATNF catalogue available in http: //www. atnf . csiro. au/research/pulsar/psrcat/ 

(Manchester et al. 2005). 



PSR 

Name 

Per 

sec 

p 

AGE 
(xlO® Yr) 

BSURF 

(Xl0l2 Q) 

E 

(xlO^^ ergs/s) 

1 

B0301+19 

1.387584 

1.30 

17.0 

1.36 

1.91 

2 

B0525+21 

3.745539 

40.1 

1.48 

12.4 

3.01 

3 

J0546+2441 

2.843850 

7.65 

5.89 

4.72 

1.31 

4 

B0656+14 

0.384891 

55.0 

0.111 

4.66 

3810 

5 

B0751+32 

1.442349 

1.08 

21.2 

1.26 

1.42 

6 

B0823+26 

0.530661 

1.71 

4.92 

0.964 

45.2 

7 

B0834+06 

1.273768 

6.80 

2.97 

2.98 

13.0 

8 

B0919+06 

0.430627 

13.7 

0.497 

2.46 

679 

9 

B0950+08 

0.253065 

0.23 

17.5 

0.244 

56.0 

10 

B1133+16 

1.187913 

3.73 

5.04 

2.13 

8.79 

11 

B1237+25 

1.382449 

0.96 

22.8 

1.17 

1.43 

12 

J1740+1000 

0.154087 

21.5 

0.114 

1.84 

23200 

13 

B1737+13 

0.803050 

1.45 

8.77 

1.09 

11.1 

14 

J1910+0714 

2.712423 

6.12 

7.02 

4.12 

1.21 

15 

B1910+20 

2.232969 

10.2 

3.48 

4.82 

3.61 

16 

B1919+21 

1.337302 

1.35 

15.7 

1.36 

2.23 

17 

B1929+10 

0.226518 

1.16 

3.10 

0.518 

393 

18 

B1944+17 

0.440618 

0.0241 

290 

0.104 

1.11 

19 

B2002+31 

2.111265 

74.6 

0.449 

12.7 

31.3 

20 

B2016+28 

0.557953 

0.148 

59.7 

0.291 

3.37 

21 

B2020+28 

0.343402 

1.89 

2.87 

0.816 

185 

22 

B2034+19 

2.074377 

2.04 

16.1 

2.08 

0.902 

23 

B2110+27 

1.202852 

2.62 

7.27 

1.80 

5.95 

24 

B2315+21 

1.444653 

1.05 

21.9 

1.24 

1.37 
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Table 3. The median values for microstructures for all the Stokes parameters. In column 2 
the longitude range is the window that has been used for analysis and this region of the pulse can 
be seen as gray regions in the average profiles given in the appendix. In column 5, 6 and 7 the 
median value is given and the values in the parenthesis correspond to the median absolute 
deviation of the P^ values. Column 8 gives the value of the bandwidth h used for analysis and 
column 9 gives the number of points in the pulse window (see text for further details). 


Median P^ from ACT histogram 


PSR 

Longitude 

Preq. 

Period 

P^ 

/i 

pL 

pV 

/i 


NPT 

Name 

Range (°) 

(GHz) 

(sec) 

(msec) 

(msec) 

(msec) 

h 


B0301+19 

-2.2 - 

-2.5 

1.2 

1.387 

1.46 (0.57) 

1.49 (0.57) 

1.25 (0.35) 

0.10 

301 

B0301+19 



1.2 

1.387 

1.19 (0.44) 

1.25 (0.48) 

1.13 (0.35) 

0.075 

301 

B0301+19 



1.2 

1.387 

1.07 (0.27) 

1.01 (0.35) 

1.07 (0.35) 

0.05 

301 

B0301+19 

-11.8 

— -5.2 

1.2 

1.387 

1.73 (0.97) 

1.85 (0.88) 

2.44 (1.10) 

0.10 

427 

B0301+19 



1.2 

1.387 

1.49 (0.62) 

1.64 (0.71) 

1.96 (0.79) 

0.075 

427 

B0301+19 



1.2 

1.387 

1.22 (0.40) 

1.31 (0.40) 

1.67 (0.66) 

0.05 

427 

B0525+21 

-2.85 

— 2.86 

1.2 

3.745 

4.17 (2.03) 

4.94 (2.47) 

6.01 (2.74) 

0.10 

1001 

B0525+21 



1.2 

3.745 

3.57 (1.76) 

4.23 (1.85) 

4.70 (2.21) 

0.075 

1001 

B0525+21 



1.2 

3.745 

2.98 (1.24) 

3.21 (1.50) 

3.51 (1.68) 

0.050 

1001 

B0525+21 

-14.3 

—-10.9 

1.2 

3.745 

5.65 (3.88) 

6.49 (3.44) 

7.50 (3.27) 

0.10 

1201 

B0525+21 



1.2 

3.745 

4.52 (2.47) 

5.48 (2.91) 

6.19 (3.27) 

0.075 

1201 

B0525+21 



1.2 

3.745 

3.39 (1.68) 

3.99 (2.21) 

4.46 (2.74) 

0.050 

1201 

J0546+2441 

-2.13 

— 1.81 

1.2 

2.843 

2.74 (1.2) 

3.04 (1.60) 

4.23 (1.85) 

0.10 

504 

J0546+2441 



1.2 

2.843 

2.38 (0.71) 

2.44 (1.24) 

3.63 (1.85) 

0.075 

504 

J0546+2441 



1.2 

2.843 

2.14 (0.88) 

1.99 (0.84) 

2.68 (1.15) 

0.050 

504 

B0656+14 

-0.61 

— 4.95 

1.2 

0.384 

0.42 (0.18) 

0.48 (0.18) 

0.72 (0.27) 

0.10 

101 

B0656+14 



1.2 

0.384 

0.42 (0.09) 

0.42 (0.09) 

0.66 (0.27) 

0.075 

101 

B0656+14 



1.2 

0.384 

0.36 (0.09) 

0.36 (0.09) 

0.60 (0.27) 

0.05 

101 

B0751+32 

-2.24 

— 2.43 

1.2 

1.442 

2.58 (1.33) 

3.12 (1.33) 

3.66 (1.60) 

0.10 

304 

B0751+32 



1.2 

1.442 

2.10 (0.98) 

2.64 (1.33) 

3.42 (1.65) 

0.075 

304 

B0751+32 



1.2 

1.442 

1.86 (0.89) 

1.86 (1.20) 

2.25 (1.33) 

0.050 

304 

B0823+26 

-1.65 

— 1.41 

1.2 

0.530 

0.84 (0.53) 

0.78 (0.53) 

0.90 (0.45) 

0.10 

264 

B0823+26 



1.2 

0.530 

0.72 (0.45) 

0.72 (0.36) 

0.84 (0.45) 

0.075 

264 

B0823+26 



1.2 

0.530 

0.60 (0.27) 

0.60 (0.27) 

0.72 (0.27) 

0.050 

264 





Table 3—Continued 


Median from ACF histogram 


PSR 

Longitude 

Freq. 

Period 

P^ 

^ /i 

pL 

py 


NPT 

Name 

Range (°) 

(GHz) 

(sec) 

(msec) 

(msec) 

(msec) 

h 


B0834+06 

-1.32 — 1.46 

1.2 

1.275 

0.78 (0.35) 

0.90 (0.46) 

1.14 (0.53) 

0.10 

167 

B0834+06 


1.2 

1.275 

0.72 (0.18) 

0.84 (0.36) 

1.02 (0.45) 

0.075 

167 

B0834+06 


1.2 

1.275 

0.60 (0.18) 

0.66 (0.27) 

0.87 (0.49) 

0.050 

167 

B0834+06 

3.7 — 7.2 

1.2 

1.275 

1.14 (0.62) 

1.26 (0.71) 

2.22 (1.51) 

0.10 

210 

B0834+06 


1.2 

1.275 

0.96 (0.45) 

1.02 (0.53) 

1.44 (0.80) 

0.075 

210 

B0834+06 


1.2 

1.275 

0.90 (0.45) 

0.84 (0.36) 

1.17 (0.80) 

0.050 

210 

B0919+06 

-7.5 — 4.74 

1.2 

0.430 

0.96 (0.27) 

0.90 (0.27) 

1.83 (0.93) 

0.10 

249 

B0919+06 


1.2 

0.430 

0.90 (0.27) 

0.84 (0.27) 

1.44 (0.71) 

0.075 

249 

B0919+06 


1.2 

0.430 

0.78 (0.22) 

0.78 (0.27) 

1.02 (0.45) 

0.050 

249 

B0950+08 

-12.77 — 10.84 

1.2 

0.253 

0.66 (0.31) 

0.66 (0.31) 

0.72 (0.27) 

0.10 

279 

B0950+08 


1.2 

0.253 

0.63 (0.22) 

0.66 (0.27) 

0.75 (0.31) 

0.075 

279 

B0950+08 


1.2 

0.253 

0.54 (0.18) 

0.60 (0.18) 

0.66 (0.27) 

0.050 

279 

B1133+16 

-1.47 — 3.39 

1.2 

1.187 

1.01 (0.62) 

0.95 (0.44) 

1.13 (0.53) 

0.10 

271 

B1133+16 


1.2 

1.187 

0.83 (0.44) 

0.83 (0.35) 

1.01 (0.44) 

0.075 

271 

B1133+16 


1.2 

1.187 

0.71 (0.27) 

0.71 (0.27) 

0.83 (0.35) 

0.050 

271 

B1133+16 

3.39 — 8.8 

1.2 

1.187 

1.13 (0.71) 

1.13 (0.80) 

1.19 (0.71) 

0.10 

301 

B1133+16 


1.2 

1.187 

0.95 (0.62) 

0.95 (0.53) 

1.01 (0.53) 

0.075 

301 

B1133+16 


1.2 

1.187 

0.66 (0.27) 

0.71 (0.27) 

0.83 (0.44) 

0.050 

301 

B1237+25 

-2.07 — 4.3 

1.2 

1.382 

3.06 (1.60) 

3.30 (1.42) 

3.12 (1.07) 

0.01 

401 

B1237+25 


1.2 

1.382 

2.16 (1.16) 

2.34 (1.16) 

2.58 (1.07) 

0.075 

401 

B1237+25 


1.2 

1.382 

1.38 (0.53) 

1.44 (0.71) 

2.10 (1.07) 

0.050 

401 

B1237+25 

8.77 — 11.87 

1.2 

1.382 

1.38 (0.45) 

1.50 (0.71) 

1.98 (0.98) 

0.01 

201 

B1237+25 


1.2 

1.382 

1.26 (0.45) 

1.26 (0.45) 

1.98 (0.98) 

0.075 

201 



Table 3—Continued 


Median from ACF histogram 


PSR 

Longitude 

Preq. 

Period 

P^ 

pL 

pV 


NPT 

Name 

Range (°) 

(GHz) 

(sec) 

(msec) 

(msec) 

(msec) 

h 


B1237+25 



1.2 

1.382 

1.20 (0.53) 

1.20 (0.53) 

1.53 (0.76) 

0.050 

201 

J1740+1000 

-4.87 

— 10.7 

1.2 

0.154 

0.54 (0.18) 

0.54 (0.18) 

0.93 (0.44) 

0.10 

113 

J1740+1000 



1.2 

0.154 

0.48 (0.18) 

0.48 (0.18) 

0.78 (0.36) 

0.075 

113 

J1740+1000 



1.2 

0.154 

0.42 (0.09) 

0.48 (0.09) 

0.60 (0.18) 

0.050 

113 

B1737+13 

-7.94 

— 5.41 

1.2 

0.803 

1.83 (1.07) 

1.95 (0.76) 

2.52 (1.20) 

0.10 

501 

B1737+13 



1.2 

0.803 

1.80 (0.80) 

1.71 (0.76) 

2.13 (0.89) 

0.075 

501 

B1737+13 



1.2 

0.803 

1.44 (0.58) 

1.47 (0.58) 

1.92 (0.71) 

0.050 

501 

J1910+0714 

-1.41 

— 3.91 

1.2 

2.712 

5.71 (1.85) 

5.65 (3.18) 

6.37 (2.82) 

0.10 

675 

J1910+0714 



1.2 

2.712 

5.36 (2.21) 

4.64 (2.38) 

5.65 (3.27) 

0.075 

675 

J1910+0714 



1.2 

2.712 

4.40 (2.56) 

4.23 (2.29) 

3.81 (1.59) 

0.050 

675 

B1910+20 

-1.21 

— 6.46 

1.2 

2.233 

5.74 (3.35) 

5.06 (2.12) 

5.65 (2.03) 

0.10 

801 

B1910+20 



1.2 

2.233 

4.40 (2.12) 

4.40 (1.85) 

4.82 (1.76) 

0.075 

801 

B1910+20 



1.2 

2.233 

3.78 (1.99) 

3.78 (1.81) 

4.37 (2.03) 

0.050 

801 

B1919+21 

-2.45 

— 3.10 

1.5 

1.337 

1.43 (0.71) 

2.83 (1.37) 

2.92 (1.50) 

0.10 

348 

B1919+21 



1.5 

1.337 

1.31 (0.62) 

2.32 (1.15) 

2.47 (0.97) 

0.075 

348 

B1919+21 



1.5 

1.337 

1.13 (0.35) 

1.61 (0.71) 

1.93 (0.93) 

0.050 

348 

B1919+21 

5.28 - 

- 8.48 

1.5 

1.337 

2.20 (1.50) 

2.23 (0.84) 

2.62 (1.15) 

0.10 

201 

B1919+21 



1.5 

1.337 

1.55 (1.15) 

1.96 (0.88) 

1.96 (0.88) 

0.075 

201 

B1919+21 



1.5 

1.337 

0.77 (0.26) 

1.85 (0.26) 

1.01 (0.53) 

0.050 

201 

B1929+10 

-11.1 

— 7.84 

1.2 

0.226 

0.54 (0.26) 

0.54 (0.26) 

0.54 (0.26) 

0.10 

201 

B1929+10 



1.2 

0.226 

0.48 (0.18) 

0.48 (0.18) 

0.48 (0.18) 

0.075 

201 

B1929+10 



1.2 

0.226 

0.42 (0.09) 

0.42 (0.18) 

0.48 (0.18) 

0.050 

201 

B1944+17 

-9.67 

— 4.91 

1.2 

0.440 

0.77 (0.27) 

0.77 (0.35) 

1.13 (0.57) 

0.10 

301 



Table 3—Continued 


Median from ACF histogram 


PSR 

Longitude 

Preq. 

Period 

P^ 

pL 

pi/ 


NPT 

Name 

Range (°) 

(GHz) 

(sec) 

(msec) 

(msec) 

(msec) 

h 


B1944+17 



1.2 

0.440 

0.71 (0.18) 

0.77 (0.27) 

0.95 (0.44) 

0.075 

301 

B1944+17 



1.2 

0.440 

0.71 (0.18) 

0.71 (0.18) 

0.77 (0.27) 

0.050 

301 

B2002+31 

-2.48 - 

- 2.01 

1.2 

2.111 

1.31 (0.71) 

1.73 (0.71) 

2.20 (0.62) 

0.10 

445 

B2002+31 



1.2 

2.111 

1.25 (0.53) 

1.67 (0.62) 

2.02 (0.62) 

0.075 

445 

B2002+31 



1.2 

2.111 

1.01 (0.35) 

1.43 (0.44) 

1.85 (0.62) 

0.050 

445 

B2016+28 

-5.34 - 

- 4.45 

1.2 

0.557 

0.66 (0.18) 

0.66 (0.18) 

0.86 (0.40) 

0.10 

256 

B2016+28 



1.2 

0.557 

0.66 (0.18) 

0.66 (0.18) 

0.83 (0.27) 

0.075 

256 

B2016+28 



1.2 

0.557 

0.66 (0.18) 

0.66 (0.18) 

0.71 (0.27) 

0.050 

256 

B2020+28 

-11.27 

— -5.03 

1.2 

0.343 

0.66 (0.27) 

0.66 (0.27) 

0.60 (0.18) 

0.10 

101 

B2020+28 



1.2 

0.343 

0.60 (0.27) 

0.60 (0.27) 

0.54 (0.18) 

0.075 

101 

B2020+28 



1.2 

0.343 

0.48 (0.18) 

0.48 (0.18) 

0.48 (0.18) 

0.050 

101 

B2020+28 

-3.78 - 

- 1.21 

1.2 

0.343 

0.54 (0.27) 

0.60 (0.27) 

0.60 (0.27) 

0.10 

98 

B2020+28 



1.2 

0.343 

0.48 (0.18) 

0.54 (0.18) 

0.54 (0.18) 

0.075 

98 

B2020+28 



1.2 

0.343 

0.42 (0.09) 

0.42 (0.18) 

0.48 (0.18) 

0.050 

98 

B2034+19 

-5.25 - 

--1.49 

1.2 

2.074 

3.36 (1.24) 

3.12 (1.10) 

3.04 (1.24) 

0.10 

361 

B2034+19 



1.2 

2.074 

3.10 (1.50) 

2.86 (1.15) 

2.92 (1.37) 

0.075 

361 

B2034+19 



1.2 

2.074 

2.80 (1.68) 

2.62 (1.06) 

2.80 (1.41) 

0.050 

361 

B2110+27 

-1.3 — 

2.26 

1.47 

1.202 

1.19 (0.09) 

1.25 (0.09) 

1.19 (0.44) 

0.10 

201 

B2110+27 



1.47 

1.202 

1.13 (0.09) 

1.19 (0.09) 

1.07 (0.35) 

0.075 

201 

B2110+27 



1.47 

1.202 

1.07 (0.18) 

0.95 (0.27) 

0.83 (0.18) 

0.050 

201 

B2315+21 

-1.79 - 

- 1.17 

1.47 

1.444 

1.07 (0.09) 

1.01 (0.09) 

1.19 (0.35) 

0.10 

201 

B2315+21 



1.47 

1.444 

1.01 (0.09) 

1.01 (0.09) 

1.07 (0.27) 

0.075 

201 

B2315+21 



1.47 

1.444 

0.95 (0.09) 

0.95 (0.18) 

0.83 (0.18) 

0.050 

201 
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A. Methods 

A.l. Models for data and noise 

The measured signal Yi at time ti is modeled as 

Y^ = fiU) + (Al) 

where / embodies the true but unknown variation of the subpnlse signal with time, and represents 
the noise in the data. The N time points ti where signal is measured are assumed to be equidistant; 
be., of the form ti = i5 where i = 0,1,... , A — 1. ((^ = 59.5/iS is the time resolution for all the 
data sets analyzed in this work). Without loss of generality, we take 5 = 1 (which amounts to 
redefining the time unit). The noise variables e^ are assumed to be independent with mean 0 and 
a (constant) variance cr^. Our present analysis does not require an estimate of the noise variance 
explicitly (except for the purpose of choosing high-S/N pulses; see Sec. 3). We have modeled the 
true signal / underlying the subpulse data as 

f{t) = e{t) +m(t), (A2) 

where e and m respectively represent a smooth slow-varying envelope and faster-scale qnasiperiodic 
modulations corresponding to microstructure. The methodology in Sec. A.3 is used to estimate / 
optimally, whereas e is estimated using the method in Sec. A.4. The difference f{t) — e(t) is taken 
to be the estimated microstrncture. Fig. 4 shows examples of such fit-envelope-microstrncture 
decomposition for a few snbpulses. The model Al is inspired by the observation based on Fig. 4 
that the ACF is dominated by the envelope feature which should be removed from the data to 
extract the microstrncture. 


A.2. Outliers 

Radio-frequency interference (RFI) over timescales that are much smaller than the pulse period 
may introduce isolated outliers in the subpulse time series which may have detrimental effect on the 
downstream analysis. For example, outliers may cause fiuctnations in the fit that may have an effect 
on the ACF used for timescale identification. Therefore, before engaging into any further analysis, 
we first assess the extent of contamination in a data set (through visual inspection of snbpulses) 
and, if deemed necessary, we identify (and “curate”) potential outliers in all the subpnlse time 
series in the data set heuristically as explained below. In this work, we have used this approach 
because of its simplicity. However, we note that a better alternative to this heuristic approach 
could be the use of robust nonparametric regression methods (Hardle 1992; Efromovich 1999) that 
deal with outliers in the regression context in a unified manner, or outlier identification methods 
with well-established formal properties. 

Fig. Al shows examples of snbpulses (#42, #44, #45) with clear outliers from a PSR B0525+21 
data set. In fact, this was the only data set in the present study which showed clear outliers. 
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Further, for this data set, it was known that these outliers were caused by radio communication 
from a satellite in the field of view. Fig. Al also shows a subpulse (#298) from the same data 
set which illustrates a failure of our heuristic method for outlier identification. This subpulse has 
strong and highly ringing microstructure which results into too many data points in the ringing 
part being wrongly tagged as potential outliers. Since we expect at most a few isolated outliers 
in any subpulse (say, less than 1% of the subpulse length), such over-identification can be easily 
detected. Outlier correction, as described below, is not applied to such subpulses. 


Identifying potential outliers. When the data contain an unknown trend /, it is necessary to 
identify a potential outlier with respect to data values in its local neighborhood. Therefore, we 
first obtain a pilot fit fo through kernel regression (Wasserman 2006; Hardle 1992). The smoothing 
bandwidth for this pilot fit may be taken as a constant b times the optimal smoothing bandwidth 
obtained via cross-validation; slight oversmoothing {b > 1) is expected to help make the fit fo 
somewhat robust against outliers. Any data point i with a relatively large residual ri = Yi — foifi) 
is a potential outlier: Specifically, following the standard prescription (Frigge et al. 1989), the ith 
data point is tagged as a potential outlier if > rrs + a x (rrs — r 25 ) or Vi < r 25 — a x (rrs — r 25 ), 
where and r 25 are respectively the 75 and 25 percentile points for the residual vector — r 25 
is therefore the interquartile range). The value of the multiplier a is taken to be 3 instead of the 
conventional 1.5 to make outlier identification more conservative. 


Imputing potential outliers back into data. We replace a potential isolated outlier value Yi 
by the median of signal values in a small local neighborhood of In this work, we have taken the 
size of the local neighborhood to be 5 points on either side of ti. That is, we replace a potential 
outlier value Yi with the median of Yi_ 4 , Yi_s, Yi_ 2 , T^-i, T^+i, T^+ 2 , ^i+3, ^i+4, ^i+5, with data 
boundaries handled by truncation. 


A.3. Model-independent fit to subpulse time series data 

The immense variability in the pulse shapes suggests that the separation of probable signal 
from noise is best done using a nonparametric regression method that only makes mild assumptions 
about the true but unknown regression function without pre-specifying its mathematical form. In 
particular, we use smoothing spline regression (Green & Silverman 1994) which may be seen as a 
penalized least-squares method with a penalty on the roughness of the fit. Optimal value for the 
smoothing parameter is chosen using the standard method of (generalized) cross-validation. The 
smoothness of the resultant fit can be expressed in terms of degrees of freedom (DF). DF can be 
interpreted as the effective number of free parameters of the fit. Specifically, a fit with greater 
fluctuation has larger DF, while a smoother fit has smaller DF. 
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Pulse 45 



Pulse 298 



Fig. Al.— Subpulses 42, 44, 45, and 298 from a PSR B0525+21 data set illustrating the behavior 
of our heuristic outlier identification method. See Sec. A.2 for description. 
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A.4. Estimating the pulse envelope 

We use the Nadaraya-Watson kernel smoother (Wasserman 2006; Hardle 1992) to estimate the 
envelope of a subpulse by judiciously oversmoothing the data. Through trial-and-error, we found 
that the smoothing bandwidth H ^ 0.075 x N is more or less reasonable for all the pulsar data sets 
we analyzed. However, considering that the smoothing bandwidth for the envelope is perhaps the 
most critical parameter of our methodology, and that there does not appear to be any principled 
way of arriving at a prescription for it, we report results for a number of values of h = H/N over 
a plausible range; the rational for this range is discussed in Sec. 4. 

B. Plots for average profiles 

The average profiles and polarization position angle (PPA) histograms for all the pulsars used 
for microstructure analysis given as boldface in Table. 1 are presented. 
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B0301+19 Arecibo 1220.67 MHz 1.387 s RM: -8.300 rad m"^ 



-10 -505 

Longitude [ ° ] 


B0525+21 Arecibo 1220.67 MHz 3.745 s RM: -39.600 rod m"^ 



-15 -10 -5 0 

Longitude [ ° ] 


J0546+2441 Arecibo 1220.54 MHz 2.844 s RM: -19.200 rod m"^ 


4 = 



-6 -4 -2 0 2 4 

Longitude [ ° ] 

B0656+14 Arecibo 1220.67 MHz 0.385 s RM: 23.500 rod m“^ 



Fig. A2.— PPA histograms as in Fig. A2 for pulsars B0301+19, J0546+2441, B0525+21 and 
B0656+14 observed with time resolutions of 59.5//sec. where the instrument and band is indicated 
above each plot. The respective upper panels give the total power (black), total linear (red) and 
circular polarization LH-RH (blue). The lower panels give the polarization-angle density as black 
dots and the average PPA is plotted as red points. 























- 41 - 


B0751+32 Arecibo 1220.67 MHz 1.442 s RM: -7.000 rad m“ 


< o 


10 

Longitude [ ° ] 



B0823+26 


Arecibo 1220.67 MHz 0.531 s RM: 


5.900 rod m' 



o 

tn 


o 


< o 


Longitude [ ° ] 




Longitude [ ° ] 


Longitude [ ° ] 


Fig. A3.— PPA histograms as in Fig. A2 for pulsars B0751+32, B0823+26, B0834+06 and 
B0919+06 observed with time resolutions of 59.5yasec. 
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237+25 Arecibo 1220.67 MHz 1.382 s RM: -0.330 rod m" 


J1740+1000 Arecibo 1220.67 MHz 0.154 s RM: 23.800 rod m"^ 


o 

o 

o 

+ 


1—1 o 

o 




0 5 10 15 

Longitude [ ° ] 


Longitude [ ° ] 


Fig. A4.— PPA histograms as in Fig. A2 for pulsars B0950+08, B1133+16, B1237+25 and 
J1740+1000 observed with time resolutions of 59.5/isec. 
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J1910+0714 Arecibo 


1220.75 MHz 2.713 s 


RM: 182.000 rad m 



Longitude [ ° ] 


B1919+21 Arecibo 1563.01 MHz 1.337 s RM: -16.500 rad m 



Longitude [ ° ] 


Longitude [ ° ] 


Fig. A5.— PPA histograms as in Fig. A2 for pulsars B1737+13, J1910+0714, B1910+20 and 
B1919+21 observed with time resolutions of 59.5yasec. 
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Longitude [ ° ] 


Longitude [ ° ] 


Fig. A6.— PPA histograms as in Fig. A2 for pulsars B1929+10, B1944+17, B2002+31 and 
B2016+28 observed with time resolutions of 59.5yasec. 
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B2020+28 


Arecibo 1220.67 MHz 0.343 s RM: -74.700 rad m“^ 


B2034+19 


Arecibo 1220.63 MHz 2.075 s 


RM: -97.000 rad m“^ 



1—1 o 


-10 -5 0 5 10 

Longitude [ ° ] 




Fig. A7.— PPA histograms as in Fig. A2 for pulsars B2020+28, B2034+19, B2110+27, and 
B2315+21 observed with time resolutions of 59.5yasec. 
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Supplemental Materials: Polarized quasiperiodic structures in pulsar radio 
emission reflect temporal modulations of non-stationary plasma flow 

C. Results of pulsar microstructure analysis 

The results of our pulsar microstructure timescale analysis pipeline (as described in Sec. 3 and 
appendix of the main paper) are available for download as a tarball with the URL 
ftp: //wm. ncra. tifr. res . in/dmitra/hires^nalysis_pdf . tar. gz 

or 

http://cms.unipune.ac.in/reports/tr-20150223/ 

Untarring this tarball creates the directory hires_analysis_pdf/. This directory contains two 
pdf files for each pulsar data set analyzed, with names of the form 

<dataset>-<envelope_smoothing_bandwidth>-pulses3. pdf 
containing pnlsewise fit, envelope, microstructure, and ACT plots (where rejected pulses are indi¬ 
cated using gray shading overlaid on plots), and 

<dataset>-<envelope_sinoothing_bandwidth>-hist3. pdf 
containing microstructure timescale histograms. For example, the file B0301+19LC-0p05-pulses3.pdf 
and B0301+19LC-0p05-hist3.pdf respectively contain pnlsewise plots and microstructure timescale 
histograms for the data set B0301+19LC when analyzed using envelope smoothing bandwidth 
h = 0.05. 

D. Ascii Profiles 

The average profiles obtained for onr microstrnctnre observations can be downloaded as ascii 
profiles from ftp: //wm. ncra. tifr. res . in/dmitra/hires^scprof . tar. gz 

or 

http://cms.unipune.ac.in/reports/tr-20150223/ 

as a tarball. Untarring the tarball creates a directory “asc_profiles” which has the ascii files 
for each pulsar. A README file describes the contents of the files. 

E. Plots for average profiles 

The average profiles and polarization position angle (PPA) histograms for the pulsars which 
were not used for microstrnctnre analysis given as non-boldface pulsars in Table. 1 are presented. 
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Table SI. Additional details on the data sets analyzed (compare with Table 2, main paper). 
Column 7 is the pulse selection threshold on the fit degrees of freedom (see Sec. 3, main paper). 


No. 

Dataset 

Longitude 
Range (°) 

9 ^ of Pulses in 
the Dataset 

NPT 

Outliers 

Curated? 

Percentile 

Cut-Off 

1. 

B0301+19 

-2.2 — 2.5 

120 

301 

N 

10 

2. 

B0301+19LC 

-11.8 — -5.2 

129 

427 

N 

10 

3. 

J0546+2441 

-2.13 — 1.81 

105 

504 

N 

10 

4. 

B0525+21 

-2.85 — 2.86 

300 

1001 

Y 

25 

5. 

B0525+21LC 

-14.3 — -10.9 

300 

1201 

Y 

25 

6. 

B0656+14 

-0.61 — 4.95 

155 

101 

N 

5 

7. 

B0751+32 

-2.24 — 2.43 

84 

304 

N 

10 

8. 

B0823+26 

-1.65 — 1.41 

226 

264 

N 

10 

9. 

B0834+06LC 

-1.32 — 1.46 

100 

167 

N 

10 

10. 

B0834+06C2 

3.7 — 7.2 

100 

210 

N 

10 

11. 

B0919+06 

-7.5 — 4.74 

104 

249 

N 

10 

12. 

B0950+08 

-12.77 — 10.84 

102 

279 

N 

10 

13. 

B1133+16LC 

-1.47 — 3.39 

210 

271 

N 

10 

14. 

B1133+16T 

3.39 — 8.8 

210 

301 

N 

25 

15. 

B1237+25 

-2.07 — 4.3 

108 

401 

N 

10 

16. 

B1237+25TC 

8.77 — 11.87 

108 

201 

N 

10 

17. 

B1737+13 

-7.94 — 5.41 

93 

501 

N 

10 

18. 

J1740+1000 

-4.87 — 10.7 

116 

113 

N 

10 

19. 

B1910+20 

-1.21 — 6.46 

341 

801 

N 

10 

20. 

J1910+0714 

-1.41 — 3.91 

102 

675 

N 

10 

21. 

B1919+21LC 

-2.45 — 3.10 

85 

348 

N 

10 

22. 

B1919+21TC 

5.28 — 8.48 

68 

201 

N 

10 

23. 

B1929+10 

-11.1 — 7.84 

529 

201 

N 

10 

24. 

B1944+17 

-9.67 — 4.91 

136 

301 

N 

10 

25. 

B2002+31 

-2.48 — 2.01 

105 

445 

N 

10 

26. 

B2016+28 

-5.34 — 4.45 

341 

256 

N 

10 

27. 

B2020+28LC 

-11.27 — -5.03 

348 

101 

N 

10 

28. 

B2020+28 

-3.78 — 1.21 

580 

98 

N 

10 

29. 

B2034+19 

-5.25 — -1.49 

134 

361 

N 

10 

30. 

B2110+27 

-1.3 — 2.26 

108 

201 

N 

10 

31. 

B2315+21 

-1.79 — 1.17 

74 

201 

N 

10 
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B0301+19 Arecibo 337.67 MHz 1.388 s RM: -8.300 rad 



Longitude [ ° ] 


J0546+24 Arecibo 351.76 MHz 2.844 s RM: -19.600 rod m"^ 



-5 


Longitude [ ° ] 




-20 -10 0 10 20 
Longitude [ ° ] 


Fig. SI.— PPA histograms for pulsars PSR B0301+19, B0525+21, J0546+2441 and B1929+10 
observed with time resolution of 59.5 /rsec, where the instrument and band is indicated above 
each plot. The respective upper panels give the total power (black), total linear (red) and circular 
polarization LH-RH (blue). The lower panels give the polarization-angle density as black dots and 
the average PPA is plotted as red points. 
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B2016+28 Arecibo 337.67 MHz 0.558 s RM: -34.600 rad m"^ 


B2044+15 Arecibo 337.67 MHz 1.138 s RM: -100.000 rad m"^ 


-5 0 5 

Longitude [ ° ] 


10 15 



Fig. S2.— PPA histograms as in Fig. SI for pulsars B1919+21, B1944+17, B2016+28 and 
B2044+15 observed with time resolutions of 59.5 yusec. 


























B2110+27 


Arecibo 337.67 MHz 1.203 s RM: -37.000 rad 


B2315+21 


Arecibo 352.81 MHz 1.445 s RM: -35.800 rod m“ 
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Fig. S3.— PPA histograms as in Fig. SI for pulsars B2110+27, B2315+21 and B2020+28 observed 
with time resolutions of 59.5 yusec. 
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-50 0 50 100 
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-50 0 50 
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J1720+2150 Arecibo 1220.54 MHz 1.616 s RM: 48.000 rod m 


-5 0 5 10 15 
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Fig. S4.— PPA histograms as in Fig. SI for pulsars J1503+2111, B1534+12, J1713+0747 and 
J1720+2150 observed with time resolutions of 59.5 /rsec. 
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Fig. S5.— PPA histograms as in Fig. SI for pulsars J1746+2245, B1855+09, B1901+10 and 
B2113+14 observed with time resolutions of 59.5 //sec. 



































